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Abstract 

We study universality in the electroproduction of vector mesons using a unified nonper- 
turbative approach which has already proved to reproduce extremely well the available 
experimental data. In this framework, after the extraction of factors that are specific of 
each vector meson, we arrive at a reduced integrated elastic cross section which is universal. 
Our calculations suggest a finite infrared behavior for the strong coupling constant. 



I. INTRODUCTION 



Experiments show a striking similarity among the production cross sections of the 
various vector mesons, pointing to the existence of a universal content in these 
processes. In the present work we study this universality property using a nonper- 
turbative procedure which describes in a unified way the cross sections for elastic 
electroproduction of the S-wave vector mesons. 

We employ a general treatment of high energy scattering based on the functional 
integral approach to QCD [HEJE] an d on the WKB method. The functional integrals 
are calculated nonperturbatively, in an extended stochastic vacuum model jU [5]. 
This nonperturbative approach has been successfully applied in many fields, from 
hadron spectroscopy to high energy scattering [HI |Sj, DVCS [9 J and structure 
functions [TU] . Recently, the method has been shown to provide a good description 
of the existing observed data on vector meson electroproduction and a comparison 
with the characteristic features of the perturbative approaches [TU [121 EES HH H3 
IT6| [T7J [18] has been extensively discussed [19] . 

In our scheme the structure of the transition amplitude is defined in terms of color 
neutral quark-antiquark states (color dipoles) that displace themselves in the exter- 
nal QCD vacuum field. The fundamental elementary dynamical effect is given by the 
loop-loop interaction: Wilson loops formed by dipoles in the photon-meson overlap 
product and in the proton (for simplicity treated in a diquark model). The tran- 
sition to observable electroproduction amplitudes of hadrons is obtained through a 
superposition of the loop-loop amplitudes with the light cone wave functions of the 
hadrons and the photon used as weights. This procedure allows the calculation of 
polarized transition amplitudes and angular distributions within a unified view of 
all processes of vector meson electroproduction. We need no external inputs, and 
the only intervening quantities are inherently calculated in the nonperturbative ap- 
proach; no new free parameters have to be introduced. A more detailed description 
of the method can be found in [6]. 

In order to look for the structure of the universal behavior that appears in the 
data we concentrate on the photon virtuality (Q 2 ) dependence of the integrated 
elastic cross sections. This is the simplest observed quantity, determined almost 
fully by the shape of the photon and vector meson wave functions, both described as 
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packets of quark-antiquark dipoles. Universality is inherent to the nonperturbative 
nature of the method here applied, and is reproduced accurately by our theoretical 
calculations. 

We first consider the loop-loop amplitude which depends on the dipole sizes and 
orientation and on the impact parameter. The amplitude for the hadronic process 
in the impact parameter space is obtained by integrating the square of the loop-loop 
amplitude over the orientation and sizes of the dipoles, with the above mentioned 
weight factors, i.e., the squares of the proton wave function and the square of the 
overlaps of photon and meson wave functions. A Fourier transform leads from 
the impact parameter picture in the transverse plane to the momentum transfer 
description. After integration over the geometric weighted distribution of loops, the 
correct Q 2 dependence and normalization of cross sections emerge [20] , which turns 
out to have form o = Ay/(1 + Q 2 /M v ) n , where My is the vector meson mass. We 
notice that this essential variable arises naturaly in our calculations, in spite of My 
not being introduced a priori in this framework. Actually the vector meson mass 
appears numerically only through the expression, given in the next section, which 
relates the electromagnetic coupling fy to the measured decay rate T v ^ e + e -. 

The numerical factors Ay that fix the normalization of the cross sections, as well 
as the power n, are determined uniquely by the model, without free parameters, 
reproducing the data with good precision. The values of n differ little from one 
vector meson to another [20]. Predicted cross sections cover a range of 5 orders of 
magnitude [19]. 

The previously known property regarding universality is that, divided by the square 
of the quark charges of the vector mesons, (in the p and u meson cases, formed with 
differents quarks u and d, this means the average effective charges), the values of 
the constans Ay are quite similar. 

We here show an important improvement of this universality property: instead of 
extracting the squared charges, we extract the squared electromagnetic couplings 
fy. Besides leading to a closer grouping of the cross section values, this criterium 
to reach universality is more fundamental, as the couplings fix the properties of the 
wave functions. The electromaqnetic coupling, experimentally accessible through 
the electromagnetic decay rates, is directly related to the value of the wave function 
at the origin. The interesting comparison of the universality in cross section values 
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obtained with extraction of squared charges and with extraction of squared couplings 
is shown in section 3. 

We then go a step further, and examine the expression that relates the couplings 
fv to the experimental electromagnetic decay rates T v ^ e + e -. A first order radiative 
correction to the decay rate [2T1 122] introduces a factor 

- Q£) a s (m f ) 

with respect to the bare decay rate. The situation is that the fundamental univer- 
sality property should be found with extraction of the pure coupling, fy\ not of 
the effective one fv, with the radiative correction factor. 

The first order radiative correction to the decay rate depends on the intervening 
quark mass, because of the running property of the strong coupling constant otsirrif). 
This quantity is known for values of the argument rrif that are larger than about 
1 GeV. Thus the bare coupling can be numerically determined for the charm and 
bottom cases. We show that, using the resulting values of bare couplings in J/ip and 
T decays, perfect matching of the corresponding reduced cross sections is obtained. 

For the light mesons we have to consider effective values of the strong coupling con- 
stant. Studies about the infrared regularization of the gluon propagator point to 
finite values for as{0). We explore the possible criterion that the complete univer- 
sality in the cross sections fixes effective values for as(rrif) , and obtain interesting 
results at the light quark masses. 

Our paper is organized as follows. In Section (|llj) we very briefly expose our approach 



to the analysis of vector meson production. In Section (III) we show that, with the 
appropriate reductions, our calculations lead to a description of these processes 
which exhibits a universality determined by properties of the overlap of photon and 



vector meson wave functions. Section (IV) concludes with a summary and discussion 



of our results. In the Appendix we present details of the vector meson light cone 
wave functions. 



II. BASIC FORMULA AND RESULTS 



For convenience we present here some basic formulae and results developed in our 
previous work on photo- and electroproduction [191 EQj [231 EH ESI |26] , where details 
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can be found. 

The amplitude for electroproduction of a vector meson V in polarization state A is 
written in our framework 

T-y* p ^vp,x(s, t; Q 2 ) = J d 2 Ridz 1 p^ v ,x(Q 2 ] z 1 , Ri) J(s, q, z x , Ri) , (2.1) 

with 

J(s,q,zi,Ri) = J rf 2 R 2 c/ 2 be- iq - b |^(R 2 )| 2 S(s,6,z 1 ,R 1 ,^ 2 = 1/2, R 2 ) . (2.2) 

Here f/> p (R 2 ) is the proton wave function and 

PrvAQ 2 ; zi, Ri) = i>v\{zi, Ri)*Vva(<2 2 ; zi, Ri) (2.3) 

represents the overlap of photon and vector meson wave functions; 
S(s, b, z\, Ri, 1/2, R 2 ) is the scattering amplitude of two dipoles with sepa- 
ration vectors Ri,R 2 , colliding with impact parameter vector b [25]; q is the 
momentum transfer 

t=(pv~ P ;f = -q 2 - m 2 p (Q 2 + Ml) js 2 + 0(^ 3 ) « -q 2 (2.4) 

In the expressions above, z is the longitudinal momentum fraction of the quark in 
the virtual photon and in the vector meson and Q 2 = —p* 2 is the photon virtuality. 

The differential cross section is given by 

da 1 T| 2 . (2.5) 



d\t\ 16irs 2 

Here, s is the center of mass energy of the 7*-proton system. 



The form of overlap written in Eq. (2.3) corresponds to SCHC (s-channel helicity 



conservation); generalizations can be made introducing a matrix in the polarization 
indices. In principle a three quark wave function should be used for the proton but 
earlier investigations have shown that a quark-diquark structure leads to fenomeno- 
logically consistent results and we adopt this model [2"5] . 

The light cone wave functions of the photon and vector meson have been previously 
discussed, and are presented in the Appendix to keep the paper self contained. Here 
we write only the scalar part of the meson wave function proposed by Brodsky and 
Lepage [27], which is used in the present calculations: 



In previous papers (25J EE], another form of meson wave function has been used : 

the Bauer- Stech-Wirbel (BSW) j28j [29] . Since it results in no essential difference 
for the purposes here considered, we concentrate on the BL wave function. 

In the expression above m/ is the quark mass and z — 1 — z . We note that 4>bl has 
two parameters: the normalization constant N ( fixed by the wave function normal- 
ization) and the transverse size parameter uj. In order to determine uj [231 12H [25] we 
turn to the Van Royen-Weisskopf formula (including color) [30] which relates the 
leptonic decay width T v _> e + e - to the wave function at the origin ip(0) 

r( ^ e+e - )= is™yi 2 (2 . 7) 

where a = 1/137.036, ey is the effective charge in the vector meson, in units of the 
elementary charge e and My is the vector meson mass. This relation is the lowest 
order static approximation for the electromagnetic leptonic decay width of a vector 
meson [3~T1. 



This condition on the wave function, which determines u, is related to fy, the 
coupling of the vector meson to the electromagnetic current, defined through 

< 0|J£,(0)|V(g,A) >= efvM v e»(q,\) (2.8) 



If we use Eq. 2.8 in the S matrix element < e + (p, s)e (pf, s/)\S\V(q, A) > and 
then calculate the decay rate of the vector meson, we obtain the relation of fy 

tO Ty^ e+e - [211 

2 3MyTy^ e+e - 
JV - T—h 



^From Eqs. 2.9 and 2.7 we then obtain 

2 3MyTy^ e+e - 12 2 2 

fv= w =W eyl#))l (2 ' 10) 

relating fy , Fy^ e + e - and if)(0). 

Mesons in different states of polarization are represented by different wave functions. 
Consequentely, N and uj are different for transverse and longitudinal states. 

The relations between fy and the wave functions are [19J, [25] 



Izzuj 2 



fv = ev-riy- / -[2lu 2 (z 2 + z 2 )(Azz) 2 + mi] exp - ^ _ " (2.11) 



47T My Jo 

in the transverse case and 



= e v — V3 16 u dz zz exp J _ 2 (2.12) 



47T 







2zzuj 



in the longitudinal case. 

Numerically, we use fy = fy = fv, with fy calculated through Eq. (2.9 ). The values 
thus obtained, together with other vector meson data, are presented in Table ITTl in 



the Appendix. There, we also show in Table III our results for N, cu, and the radius 
for each meson in each state, calculated through the wave function normalization 



equations, Eqs. (2.11) and (2.12) above (with fy from Table |H|) and the mean 
value of the radius r. For the heavy quarks we take the renormalized masses m c 
= 1.25 GeV, mj = 4.2 GeV and for the light quarks we take the effective masses 
determined 1321 as m u = m d =0.2, m s = 0.3 GeV. 



For the quark-diquark proton wave function, fixing the longitudinal momentum 
fraction at z ~ \ and making a Gaussian ansatz, we have 



1 1 

2tt S n 



exp 



-2 , 



(2S P ) 



(2.13) 



where the size parameter is chosen as S p = 0.74 fm [33 J. 



In our model the energy dependence is motivated by the two-pomeron approach of 
Donnachie and Landshoff [31]. The hard pomeron leads to (s/so) ' 42 for Ri < r c and 
the soft pomeron leads to (s/so) ' 0808 for large dipoles, so that the integration over 
Ri in Eq. (2.1) is split into two parts [2S]. We use the numerical values r c w 0.22 
fm and Sq = (20 GeV) 2 , taken from 



III. UNIVERSALITY 

As said above, we have used a nonperturbative QCD framework to calculate the 
observables of vector meson photo- and electroproduction, and we have obtained 
very good description of the data in all cases PH EDI E3J ES ESI EH] . Fig. [l]shows the 
data and the theoretical calculations of the electroproduction of p,u,<p, J/ip and T 
vector mesons. 

The quantitative predictions made in a unique way for different kinds of vector 
mesons cover 5 orders of magnitude in the cross sections [19J. The nonperturbative 
calculation leads to cross sections that plotted against Q 2 + M v are nearly parallel 
straight lines, a result confirmed by experiments. Our work has shown that overlap of 
photon and vector meson wave functions alone determines large part of the behavior 
observed in these processes. This simple structure suggests a search for universality. 
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FIG. 1: Integrated elastic cross sections in vector meson electroproduction: experimental 
data and theoretical calculations with the stochastic vector model. 




The effect of factorization in the magnitudes of the cross sections is exhibited with 
the extraction of charge factors e y = 1/2 , 1/18 , 1/9 , 4/9 , 1/9 for p, u, 0, ip, 
and T, respectively. Fig. [2] shows the theoretical results for all vector mesons after 
charge reduction, with approximately equal values for all cross sections. Since it can 
be seen from Fig. [I] that the experimental values agree very well with the theoretical 
predictions the universality to be derived from the theoretical expressions reflects 
an universality in the data. The agreement of our nonperturbative approach with 
the experimental results is a consequence also of the construction of the meson wave 
function, whose shapes are fixed by the experimental values of the electromagnetic 
leptonic decay rates of vector mesons T v ^ e + e - . These quantities fix the transverse 
size parameter u contained in the wave function, as explained above [231 12H 125] . 
The coupling of the vector meson to the electromagnetic current (/y) is related to 
Ty^e+e- as shown in Eq. (2.10). The values obtained for f v are collected in Table [TT] 
in the Appendix. 

Fig. [3] shows the theoretical calculations of the electroproduction of p , u, (ft, J/ip 
and T vector mesons divided by the values of f v . We observe a remarkable univer- 
sality, closer than obtained with extraction of factors of squared charges. The small 
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FIG. 2: Extraction of charge factors e v from the integrated elastic cross sections in vector 
meson electroproduction. 




remaining splittings may be due to flavor dependent corrections in the calculations 
of the electromagnetic decay rates. We take these into account in the discussion 
that follows. The leading QCD radiative correction can be calculated in perturba- 
tion theory |22J ESI EI| and causes a significant suppression in the value of Y 



given by the zero-order equation (|2.7|), which becomes 

r(v 

leading to 



e e 



r n 



1 - 



16 
37 



a s (m f ) 



f v = ti° )2 F(m f ) « /, 



(0)2 
V 



1 - 



16 

3^ 



«sK) 



(3.1) 



(3.2) 



The flavor dependence comes through the running strong coupling constant as, 
which depends on m/. Since as is known at the mass values of the b and c quark 



masses 



a s {m b ) = 0.23 ± 0.01 , a s (m c ) = 0.33 ± 0.01 



the corrected Eq. (3.2) can be applied directly to T and J ftp. 



Fig. [i] shows the theoretical calculations of electroproduction of J / tp and T vector 
mesons with extraction of fy^ 2 factors, with remarkable coincidence of the reduced 



FIG. 3: Extraction of fy coupling factors from the integrated elastic cross sections in 
vector meson electroproduction. 



10 r? 

y* p ->(p,<o,4>,JA|'.Y) p 
X W=90 GeV 




values of the cross sections. We thus have precise universality for T and J /if) mesons. 

The simple behavior connecting J/if) and T cross sections, that we explained by 
sheer acting on wave function construction, is a manifestation of the factorization 
that shows that the whole Q 2 dependence of the cross sections is determined by the 
overlap of the photon and vector meson wave functions. In Fig. [5] we exhibit the 
accuracy of the calculation of the J /if) cross section using only the overlap strength 
of photon and J / if) wave functions (the overlap strength is formed by the integration 
over the internal variables of the qq pairs of the overlap function multiplied by r 2 ) 
|20| . The accurate calculation of the cross section with the stochastic vacuum model 
in the LHS is compared to the calculation with the overlap strengths in the RHS. 
A numerical factor 24000 takes into account the interaction of the proton with the 
QCD field, which has no Q 2 dependence. Actually the Q 2 dependence of cross 
sections is of the form 

Ay 

° (Q 2 + M 2 Y 

which is contained in the shape of the integration of products of the modified 
Bessel functions Kq and K\ with powers and Gaussians that appear in the over- 
lap strengths. 
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FIG. 4: Extraction of the bare f v coupling factors from the integrated elastic cross 
sections of J/ip and T mesons . 




FIG. 5: Data and calculation of J/tp cross section as a function of photon virtuality . 
In the RHS we show the Q 2 dependence of the overlap strength determined only by wave 
functions. 
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Good results for cross sections calculated only with these quantities are obtained 
whenever the range of the overlap region is small compared to the typical range of 
the nonperturbative interaction governing the process |20j . 

Let us next consider the light vector mesons. Here we may have some complications, 
since it may be necessary to include relativistic and higher-order radiative corrections 
to fv- Moreover, we will need the values of as for small m/, which are not yet 
determined experimentally. 

While there is a general agreement as to the high energy behavior of the strong cou- 
pling constant [381 ESI HQ] , the infrared region, where as(Q 2 ) develops singularities, 
is not determined. Perturbative calculations point to a singular behaviour , but the 
possibility of finite values at small Q 2 has been extensively studied in recent years. 
There are analytical, phenomenological and lattice calculations exploring this possi- 
bility, and many arguments and techniques have been proposed to solve the problem 
of the elimination of the infrared singularities [391 HI]- 

A phenomenological proposal to determine the infrared strong coupling constant is 
that of the so called physical coupling or effective charge, which consists in defining 
cts, e tt{Q 2 ) from a given physical observable [12]. This was analysed in connection 
to the renormalization-scheme dependence problem in perturbative field theories 
[4*3"| |4*4"]. This effective charge is a resummation in perturbation theory [331 HE1 HZ]- 
Brodsky et al. argue that a perturbatively calculable physical quantity can be used 
to define an effective charge by incorporating the entire radiative correction in its 
definition |45|. 



In line with this phenomenological proposal, for the light vector mesons we define an 
effective strong running coupling constant that absorbs the higher order corrections 



to fv and so, instead of the corrected Eq. (3.1), we write 



(3.3) 



Now the reduction factor in Eq. (]3.2|) is written as 

16 



F(m f ) 



3tt 



as,es(m f ) 



(3.4) 



and is determined by fitting that leads to universality of cross sections. The results 
obtained are shown in Table IIHI With this choice for the factors that account for 
the radiative corrections to the electromagnetic decays of the light vector mesons, 
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TABLE I: Reduction factors of T for (j), p and uj 



Meson mj(Lrev) 


F(fitted) 


(/> 0.3 


0.32 ± 0.01 


p 0.2 


0.32 ±0.01 


a; 0.2 


0.29 ±0.01 



FIG. 6: Extraction of renormalized f v coupling factors from the integrated elastic cross 
sections of electroproduction of all vector mesons: theoretical results and experimental 
data converted with same factors for each vector meson, as described in the text. 




we have the renormalized couplings f v for all mesons and then we may put the 
electroproduction cross section of all five vector mesons in a single line, as shown in 
Fig. [6j Through the fitted values of F(rrif) we obtain values of the effective strong 
coupling constant at low masses 

asM m f) ~ °- 4 > for m f ~ (°- 2 - °- 3 ) GeV - 
This result is similar to values found in the literature. 

On the theoretical side, there are many calculations leading to a region below 
as(Q 2 ) ~ 1 for low values of Q 2 . Let us explicitly mention a few of these. 
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In a classical paper [48J, studying gauge invariant Schwinger-Dyson equations, Corn- 
wall advocated the dynamical generation of an effective gluon mass, resulting in a 
gluon propagator that reaches a finite, non vanishing value in the deep infrared and 
in cts(Q 2 ) which is asymptotically free in the ultraviolet and attains a finite infrared 
value. This general picture is confirmed by using large volume lattices, for SU(2) 
[4*9] and SU(3) [50] pure Yang-Mills. Other recent studies are in agreement with 
Cornwall results |HJ ED E2J EH] • 

According to the original Cornwall proposal, we have 

as(M Zo ) _ a s (M Zo ) ft Q 2 + 4m 2 
a s m + 2vr M Zo ■ (3 - 5) 

Here, M Zo is the mass of Z used as reference, /?o = 11 — |n/, with rif — number of 

flavors, and m g is the gluon mass, which can be neglected in face of M Zo . 

This infrared finite strong coupling constant has been used to study the asymptotic 
pion form factor [51], resulting in a better agreement with the data. 

The values of the strong coupling constants that appear in the radiative corrections 
to the electromagnetic decays, which modify the effective couplings, are exhibited in 
Fig. [7j In the fi gure we compare the values here obtained with measurements and 
plot the perturbative and nonperturbative [IB] calculations. In the nonperturbative 
case we use m g = 0.6 GeV . The figure shows that the value of as, e s ~ 0.4 for 
low masses, obtained by our fitting, is in the range of values predicted by Cornwall 
|48J and by more recent work. In a later paper [55], and using a different technique 
(the functional Schrodinger equation), Cornwall estimated as(M 2 ) « 0.4 for M = 
0.6 GeV, which is again in agreement with our result. The same is true for ref. 
[52J, where as is found to be a bit below 0.4 for Q 2 ps 0.1 GeV 2 . In a recent paper 
[53], calculations of QCD effective charges performed within two different theoretical 
frameworks are shown to acquire a common value «s(0) = 0.7 ± 0.3. 

On the phenomenological side, several different analyses of scattering data sensitive 
to low momentum effects point to values of as{Q 2 = 0) in the range [4"T] . 

a s (0) « 0.7 ±0.3 

which is in agreement with our values. 

In addition, we mention that the underlying stochastic vacuum model also predicts 
finite strong coupling as ~ 0.8 in the infrared region [56J. 
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FIG. 7: Strong running coupling constant. The dashed lines are the usual perturbative 
expresions, and the solid line represents Eq. (3.5) with m g = 0.6 GeV . 




IV. CONCLUSIONS 

Our treatment of elastic vector meson electroproduction stresses nonperturbative 
aspects of the QCD calculation, showing that large part of the description of this 
process is supplied by the hadronic structure, namely, the wave functions of pho- 
tons and vector mesons and a simple diquark model for the proton. Details of 
the structure of the proton, that participates in the processes like an overall ob- 
ject, determined by its size only, do not enter in the description of the integrated 
observables. 

A full description of the framework of the calculations can be found in previous 
papers [El EDI ESI 123 EE]- The calculations, without free parameters, show that 
quite simple features are able to explain the data, covering a wide range of 5 orders 
of magnitude in the cross sections[19j. 

Each of the different mesons enter in the calculation characterized only by the masses 
and charges of its quark contents, and with their normalized wave function individ- 
ualized only by the corresponding electromagnetic decay rate (related to the value 
of the wave function at the origin). The main features of the Q 2 dependence of 
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electroproduction of vector mesons are contained in the overlap integral of the light 
cone wave functions of photons and mesons folded with the basic r 2 behaviour of 
the dipole cross sections [20] . 

In this paper we concentrate on the integrated elastic cross sections, studying its 
universality in the Q 2 dependence . The energy dependence, which is well described 
by the two-pomeron model, is not detailed here : we choose a value W=90 GeV for 
which experiments are available for all vector mesons. 

Our calculation gives a unified description of the electroproduction for all five S-wave 
vector mesons. The Q 2 dependence of the integrated elastic cross section shows a 
regular shape when plotted against the variable Q 2 + M v . The lines are almost 
parallel, and it is well known that, after extraction of factors equal to the squared 
electric charges of the constituent quarks in the vector mesons (with appropriate 
averages in the p and u cases), the reduced cross sections tend to concentrate near 
a common line. We show that this approximate universality becomes more exact 
when the extracted factor is the square of the electromagnetic coupling fy . This 
improved universality is interesting because this factor is more meaningful than the 
squared charge, since it is connected to the value of the wave function at the origin, 
entering in the determination of the basic parameter uj of the wave function, and 
fixing the size of the vector meson. 

The numerical value of the coupling fy for each vector meson is fixed by the experi- 
mental value of the electromagnetic decay rate Ty_ e + e -. This experimental value is 
affected by radiative corrections at the electromagnetic vertex, which depend on the 
masses of the intervening quarks. We remark that a more exact universality would 
be observed if these corrections were not present. 

The criterion of universality of reduced cross sections imposed for the light meson 
cases, using the same expression for the radiative correction, leads to effective values 
for the strong coupling constant at the low mass values. These values are in iteresting 
agreement with theoretical studies concerning the infrared behavior of the gluon 
propagator. 

The successful description of electroprodution cross sections, without free parame- 
ters, including all vector mesons, and the universality property that is inherent to 
our nonperturbative framework, shows that microscopic details are not important 
for the observables here studied. Only the dipole distribution in the wave functions 
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and the hadronic sizes, and the range of the interaction mediated by the nonpertur- 
bative QCD vacuum field, are relevant for the phenomenology. 
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Photon and vector meson light-cone wave functions 



APPENDIX A: PHOTON AND VECTOR MESON LIGHT-CONE WAVE 
FUNCTIONS 



The color part of the photon wave function can be treated separately and the result 
is an overall multiplicative factor y/W c . The helicity and spatial configuration part 
of the wave function ip^^Q 2 ; z,r, 6) is calculated in light-cone perturbation theory. 
In lowest order, for each polarization A, we write [131 123], [27J |571 [581 159] - 

271 

+m f 5 h , + 5- h)+ K (er) , (Al) 



ip r>1 {Q 2 ;z,r, 



iee 1 {z5 h ^ + 5- h _ - z5 h -5^ + )Ki{er) 



and 



where 



V> 7 *,-i(<3 2 ;z,r, 



'6a 



iee 1 (zS^+S* - zS^S^ + )Ki(er) 



Ztt 

+mfS h) -Sj i _Ko(er 
^3a 



zzQ 2 + m 2 . 



(A2) 



(A3) 



V , (A4) 

a = 1/137.036 and m/ is the quark mass, ej is the quark charge in units of the 
elementary charge for each flavor /; z(z) is the quark (antiquark) longitudinal mo- 
mentum fraction; h(h) is the quark (antiquark) helicity; K and K\ are the modified 
Bessel functions. In the longitudinal case there is an additional 5 2 (r) which ca be 
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dropped because the color interaction vanishes at zero transverse distance. The 
quark masses we use are m u = = 0.2, m s = 0.3, m c = 1.25 and m& = 4.2 (values 
in GeV). 

Meson wave functions are more model dependent than photon wave functions. We 
take for the vector meson the spin structure from the vector current leading to 
expressions similar to those of of the photon j23j 1571 ? ] . 

ip Vt+ i(z,r,9) = (-ie te d r (z6 h>+ 6h_-z6 h) -.Sj lj+ ) 

+m f S h}+ 8j li+ j4) V (z,r) , 
Vy,-i(z,r,0) = ( - ie~ %e d r {z8 hi+ 8- h _ - z8 h 

+m f 8 h _8 h _)4> v (z,r) (A5) 

and 

*l>v,o(z,r) = (ujAzz5 h) _ij </> v (z,r) . (A6) 

In this paper we use the functional form for the scalar function <py(z,r) proposed 
by Brodsky and Lepage 



(t>BL{z,r) 



N 



exp 



m f (z - I) 
2zzui 2 



2„2l 



exp[— 2zzuj r 



(A7) 



It contains two parameters, iV and u, which are determined, as explained in section 



2, by the normalization condition and the leptonic decay width, V 



which 



is related to fv, the coupling of the vector meson to the electromagnetic current 
Table [IT] shows vector meson data and the values obtained for fy through eq. 2^ 



Numerically, we use f\ 



fv = fvi with fy calculated through eq. 2.9 The val- 



ues thus obtained, together with other vector meson data, are presented in Ta- 
ble [IIJ Table III shows our results for u, N, and radius r for each meson in each 



state, calculated with eqs. 2.11 and 2.12 , the wave function normalization equation 
Jq 1 dz f d 2 r \ipv,\( z i r )| 2 = 1 an d the mean value of r. 

After summation over the helicity indices, the overlaps of the photon and vector 



meson wave functions (eq. 2.3) are 



Prv,±i(Q 2 ; z,r) = ^-^(4e/ uj 2 rzz[z 2 + z 2 ]Ki(e f r) +m 2 f K (e f rfj BL (z,r(A8) 
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TABLE II: S-wave vector meson data. The coupling fy and the electromagnetic decay 
width T v ^ e + e - are related through eq 2.9 The quantity ey is the effective quark charge 



in units of the elementary charge, determined by the qq structure of each meson. 



Meson 


My(MeV) 


ey 


T v ^ e+e - (keV) 


fy (GeV) 


p(770) 


775.49 ± 0.34 


l/y/2 


7.04 ± 0.06 


0.15645 ± 0.0043 


^(782) 


782.65 ±0.12 


1/3 


0.60 ±0.02 


0.04588 ± 0.0008 


0(1020) 


1019.455 ±0.020 


-1/3 


1.27 ±0.04 


0.07619 ±0.0012 


J/if>(lS) 3096.916 ±0.011 


2/3 


5.55 ±0.141 


0.27759 ±0.0035 


T(15) 


9460.30 ± 0.26 


-1/3 


1.340 ±0.018 


0.23839 ±0.0016 



TABLE III: Parameters of the vector meson wave functions 



Transverse Longitudinal 

w(GeV) N radius (fm) w(GeV) N radius(fm) 



p(770) 


0.2809 


2.0820 


1.0637 


0.3500 


1.8366 


0.6402 


w(782) 


0.2618 


2.0470 


1.1110 


0.3088 


1.8605 


0.7194 


0(1020) 


0.3119 


1.9201 


0.8569 


0.3654 


1.9191 


0.5973 


J/^ilS) 


0.6452 


1.4752 


0.3418 


0.7140 


2.2769 


0.2898 


T(1S) 


1.3333 


1.1816 


0.1551 


1.3851 


2.7694 


0.1456 



for the transverse case, 



3ct 

Prv,o(Q 2 ; z,r) = -16e v ^- — u z 2 z 2 Q K (e f r)4> BL (z,r) 

Z7T 



(A9) 



for the longitudinal case. 

The overlap strengths, formed by integration over the internal variables of the quark- 
antiquark pairs of the overlap function multiplied by r 2 are 

(MO) 



Y r y T (Q ) = J dz j d r r p r y±i(Q ;z,r) = e v Y r y T (Q ) 
for transverse polarization and 



Y r y L {Q' 2 ) = j dz j d 2 r r 2 p r y (Q 2 ; z,r) = e v Y r y L (Q 2 ) 



(AH) 
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for 



longitudinal polarization. 
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